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Remarks Concerning the Study of Four-Jet Events 
from Hadronic Decays of the Z^. 

H. Jeremie^, P.Leblanc, E.Lefebvre, D.Davignon, G.Karapetian; Universite de Montreal. 

Abstract : The angular correlations of four-jet events from hadronic decays of the Z° 
have been studied in the past mainly to extract from them the fundamental constants 
- - - of quantum chromodynamics called colour factors. Previous studies have used all the 
Q ■ available phase-space in order to maximize statistics. In this note we want to point out the 
O ! possibility that significant differences between experiment and theory in restricted regions 



of phase-space might have escaped detection. Such differences could be a harbinger of the 
existence of new particles. Some preliminary results are presented. 
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cN 1 Introduction 



^ ■ In second-order QCD perturbation theory four-jet events arise from the production of 
■ four hard partons and their subsequent fragmentation. The cross-section for four-parton 
O I production by electron-positron annihilation can then be written as: 



where Dqqqq(x) and Dqqgg(x) are differential, normalized, theoretical distributions of 
variables x for events whose final state consists of four quarks, and events consisting of 
^ ■ two quarks and two gluons respectively ^, while da{x) four-parton/ dx is the normalized 
distribution of the data. The value for the fraction /q is a fundamental prediction of 
^ [ QCD, it is approximately proportional to Tji, one of the QCD colour factors. It can 
^ I be extracted from a sample of four-jet events by varying /q to obtain a best fit to a 
measured distribution, if Dqqqq(x) and Dqqgg(x) are sufficiently different from each other 
(see f.ex. figure 3d) Such differences arise because in the case of qqqq events the primary 
quark-antiquark pair radiates a polarized gluon which splits into two quarks of spin 1/2, 
while in the qqgg case this intermediate gluon splits into two gluons of spin one. Angular 
momentum conservation then requires the distribution of the final quark-antiquark pair 
with respect to the initial pair to be different from that where the final partons are 
gluons. For practical reasons we preferred to work directly with /q than with Tr, the two 
approaches being equivalent. 

If the intermediary gluon would have the possibility to split into fermions other than 
those with the canonical QCD flavours, the value of /q (or Tr) would show an increase, 
constituting thus a possible signal of some new physics [2]. 



^ e-mail: j eremie@lps . umontreal . ca 

^this is the English version of an article which appeared in French in ref. [1] 

'^r'qqgg(x) contains events where both quarks radiate a gluon as well as events where one radiated 
gluon splits into two gluons. 
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The main theoretical distributions we used here were obtained from the second-order 
{0{al)) QCD matrix element calculation by Ellis, Ross and Terrano (called ERT model 
hereafter) [3], as implemented in the JETSET simulation package [4], which includes 
hadronization of the four partons. The parameters of the model were adjusted so as to re- 
produce energy-energy correlations and event-shape distributions measured by OPAL [5] ^, 
and the scaled invariant mass cut-off of the ERT four-parton generator was t/min=0.01. 

Higher order QCD processes will also contribute to the production of four-jet events, in 
which case the value of /q derived from eq.(l) will no longer yield a correct measurement of 
this quantity. However, /q still provides a useful method of comparison between data and 
models. Since this method relies only on comparisons of shapes of angular correlations 
and no actual four-quark events are ever identified, we call the measured quantity the 
"apparent" fraction of four quark events. 

To gauge the influence of higher orders, a hybrid model combining matrix elements and 
parton showers, and a matrix element model generating pure five-parton events were also 
investigated (see refs.[6] and [7]). Such a hybrid model is expected to simulate a number 
of higher order effects, while the 5-parton events are one of the important contributions 
of next-to-leading order(NLO) calculations. 

Using all events produced in the entire region of phace-space (except for a threshold) , 
previous investigations (refs. [8] to [15]) found good agreement with QCD predictions for 
the overall global value of Tr, equivalent to /q. 

As previously mentioned, in this work we will try to obtain more information by ex- 
amining the evolution of /q as a function of variables which permit a subdivision of the 
phase-space. We chose the variable + m^, the sum of the masses of the two least 
energetic jets. 

2 Experimental Method 
2.1 Event Selection 

The events were recorded with the OPAL detector [16] at LEP^, details of the trigger 
selection and on-hne flltering system can be found in refs. [17] and [18] . A standard selec- 
tion of multihadronic events was applied [19]. For the present analysis charged particles, 
measured in the central tracking system of the OPAL detector, and clusters formed by 
showers in the electromagnetic and hadronic calorimeter were used. Additional correc- 
tions were applied to minimize double counting of charged particles in the central detector 
and the calorimeters [20] . 

^The values of the fragmentation parameters we used were: PARJ(21)=0.49, PARJ(41)=1.8, 
PARJ(42)=0.6. 

"''Although data taken with the OPAL detector were used, the analysis itself was done independently 
from the OPAL collaboration. Only the authors cited above assume responsibility for this analysis. 
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2.2 Reconstruction of Four Jets 

We assigned the observed particles to jets, whose directions reflect approximately the 
directions of the hard partons presumed to be emitted before hadronization^. In this 
analysis we employ the DURHAM jet-finder [21], which uses a jet resolution parameter 
between a pair of jets defined as 

Vij = 4r ■ ^^i^h - cos%) (2) 

where Ei,Ej correspond to the energies of jets i and j, while 6ij is the angle between them, 
and Eyis is the visible energy of the whole event. 

The jet-finding algorithm combines iteratively each pair of particles k and / with the 
smallest i/ki into a new pseudo-particle by adding the individual four-momenta, until 
exactly four pseudo-particles are reconstructed, which are defined to be four jets. Then 
the minimum value of yij of all the possible combinations of the four jets is determined. 
This value is called and is used to preclassify the events^. If one stops the iteration 
of the jet-finding algorithm at five jets, then the corresponding value of yij is called y^^. 
The usual way to choose a sample of four jet events is to specify a fixed value of a 
parameter called ?/cut and select events such that y^^ < ycnt < U^- The quantities 
y"^^ and y^^ represent those values of ^cut where the event makes a transition from a 
five-jet classification to a four-jet classification, and four-jet to three-jet, respectively. 
This method of selecting four-jet events is not sufficient to exclude five- and more parton 
events reconstructed as four-jet events. We therefore accepted an event if > 0.012, but 
demanded < 0.006. This eliminates a large fraction of unresolved five-parton events 
in the sample, without reducing appreciably the number of four-jet events (22% for data 
but only 8% for ERT, at the detector level ) 

From an analysis of 3.3x10^ accepted hadronic events by OPAL at the Z° peak we 
extracted 8.4x10^ four-jet events selected as described above. 

A sample of 14.3x10*^ simulated four-jet events was used to correct for detector and 
acceptance effects. They were generated with the ERT matrix element calculation we 
want to test and subjected to the same cuts as the data. 

2.3 Measurement of the Bengtsson-Zerwas correlation 

We used as angular variable mostly the Bengtsson-Zerwas [22] correlation, other correla- 
tions arc also possible [23]. It measures the angle xbz between the two vectors x 
and ^3 X jf^. Here 'p' ^ to 'jf^^ arc the momentum vectors of jets 1 to 4, where the num- 
bering refers to jets ordered with respect to their measured energy ( Ei > E2 > E^ > £'4 ) . 

^The Monte Carlo generator model generates first partons, then fragments them into particles and 
finally tracks these particles through the detector. These three stages will be referred to in the text as 
parton level, particle level and detector level respectively. 

^This definition of y^'^ is not to be confused with a jet resolution parameter 2/34 limited to two individual 
jets numbered three and four! 
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Ordering increases the probability that jets 1 and 2 do indeed originate from the primary 
quarks as required. To minimize the sensitivity to exchanges between jets 1 and 2 or 3 
and 4, we used the symmetrized version of the B-Z correlation 

Xbz = ^( P 1 X p 2) , ( P 3 X P 4) ; if Xbz > -, Xbz = TT - xbz (3) 

so that effectively xbz runs from zero to |. To obtain well defined planes, the angles 
between jets 1 and 2, and between jets 3 and 4, were required to be smaller than 160° 
The number of accepted events quoted above includes this cut. The Bengtsson-Zerwas 
correlation discriminates well between qqqq and qqgg events, since the polarization of the 
intermediate gluon [24] influences directly the distribution of its decay products, either 
two quarks or two gluons. It measures only angles and it relies on jet energies only in so 
far as energy ordering is required. The angular resolution, as obtained from comparing 
the ERT Monte Carlo at the particle level with that at the detector level, was 14 degrees 
RMS. 



3 Results 
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Figure 1: Comparison of the distribu- 
tions of 777-3 + between real data 
and ERT simulation (normalised with 
respect to each other). An overall cor- 
rection of 1 GeV has been applied to the 
scale of the ERT distribution in order to 
make the two distributions coincide bet- 
ter. The hashed histogram corresponds 
to the distribution of four-quark events 
(qqqq), while the dotted curve repre- 
sents the qqgg events. 



The variable which we will use to subdivide the phace-space is 7773 -|- 7774, the sum 
of the intrinsic masses of the two least energetic jets. Before proceeding to calculate the 
correlations and fit them to the data, we have to assure ourselves that the raw distributions 
of the events as a function of m3 + 777,4 are satisfactorily reproduced by the ERT simulation. 
Figure 1 shows that this is the case. 

In figure 2a one can find the results for fg (full circles), obtained by adjusting f^ (i.e. the 
relative proportions of the theoretical qqqq and qqgg correlations) to obtain a best fit to 
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the experimental data®. One sees that the data follow approximately the ERT predictions 
(histogram) until 10 GeV, for values above 10 GeV there are substantial differences. 

Figure 2: a.) Comparison of the distri- 
bution of /q as a function of rris + 1714 
for the data (full circles) with the ERT 
prediction (histogram). 

b. ) Comparison of the distribution of /q 
as a function of ms + for a mixture 
of 5 parton events (15%) and 4 parton 
events (85%) (full circles) with the ERT 
prediction (histogram). The proportions 
of the mixture have been adjusted so 
as to reproduce the global (integrated) 
value of i?4q = 1.32. 

c. ) Comparison of the distribution of /q 
as a function of + 777,4 for the hybrid 
model (full circles) with the ERT predic- 
tion (histogram). Here also the propor- 
tion of the hybrid events has been ad- 
justed so as to reproduce the global (in- 
tegrated) value of i?4q — 1.32. All errors 
are statistical only. 

If we integrate over all events, the global fit yields an apparent fraction of 
f^^?' = 0.108±0.009 compared to a theoretical value of = 0.082±0.001, 

(errors statistical). 

A convenient way of comparing experiment to theory is to form the ratio 
R4q_^ju^, which yields 

i?4q = 1.32±0.10 (x^/dgf=0.9) for all accepted events. The cited error being purely 
statistical, such a value of R^q represents an approximate agreement between data and 
theory. 

But if we divide the available phase-space into two regions, one with ms + m^ <10 GeV 

and the other with 7773 + 7774 > 10 GeV, we find 

for ms + m4 < 10 GeV: i?4q = 0.88±0.11 (xVdgf=0.24), while 

for 1713 + 1714 > 10 GeV: /?4q = 2.98±0.35 (x^/dgf=1.3), errors statistical only. 

These rather large differences between the values of R4q for the two regions of phase- 
space constitute the most important result of this analysis. 

^Some bins in the region of high ms -|- 1714 in figure 2a of ref.[l] have been combined, yielding the 
results displayed here 




6 



Let us now form the double quotient 
Q = g^'^ihri^lni = 3.39±0.57, 

which is yet another convenient way to represent our results. 
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+0.18 


f.) y^5 < 0.009 


±0.42 


g.) 034 < 140 


±0.15 


h.) without correction of the ERT scale 


-0.05 


i.) configuration CH+EM, without hadronic clusters 


+0.33 


k.) JADE algorithm [25] instead of DURHAM 


-1.13 


1.) N-R correlation [23] instead of B-Z 


-0.07 


m.) distribution instead of B-Z 


+0.10 


n.) constrained fit : ^Ei — Et^t '■, J2^i — 


-0.21 


0.) statistical error 


±0.57 


total systematic error 


+0.79 
-1.30 


total error 


±1.2 



Table 1: Systematical errors of the quotient Q 



An estimation of the sytematical errors of this quantity can be found in table 1. In 
this table we consider the variations of the following quantities: 

a. ) increase of the jet resolution criterium, y^^, to 0.015; 

b. ) elimination of jets which point approximately in the direction of the beam pipe; 

c. ) elimination of jets with few particles; 

d. ) and c.) elimination of small and large angles for Xbz] 

f. ) admission of a larger number of 5-jet events; 

g. ) limitation on 6*34, the angle between jet 3 et jet 4, to 140 degrees instead of 160; 

h. ) without the correction of 1 GeV for the Monte Carlo scale; 

i. ) jets calculated without using hadronic clusters and without correction for double 

counting of energy; 

k) using the old JADE algorithm (see ref. [25]) with less jet finding resolution; 
1.) replace Xbz by the Nachtmann-Reiter correlation (see ref. [23]); 
m) replace xbz by 6'34; 

n.) fit with constraint by imposing conservation of total energy and momentum. 
One finds: 

Q = 3.391? [42, so that Q becomes approximately 
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Q = 3.4±1.2. 

The difference between experiment and theory, expressed by the value of Q>1, is 
approximately an effect of two standard deviations. 

Figure 3: a.) Comparison of the dis- 
tribution of xbz between data (full cir- 
cles) and ERT prediction (histogram) for 
all accepted events (before fitting for the 
best value of fg). 




b. ) Comparison of the distribution of 
Xbz between data (full circles) and ERT 
prediction (histogram) for events with 
7713 + 1714 < 10 GeV. 

c. ) Comparison of the distribution of 
Xbz between data (full circles) and ERT 
prediction (histogram) for events with 
1713 + 1714 > 10 GeV. 

d. ) Comparison of the xbz distribu- 
tion for the four-quark events (full cir- 
cles) with the one for two quark and two 
gluon events (histogram), as calculated 
with the ERT prediction for all accepted 



c) Xbz (degrees) d.) Xbz (degrees) CVCUtS. NotC the differCUCC iu SCalc for 

this last graph compared to the three 
previous ones. 

All the distributions are normalised with 
repect to each other, they are not cor- 
rected for detector effects and all errors 
are only statistical. 

In figure 3 we display directly the xbz distributions which yielded the above-mentioned 
results. In figure 3a the experimental and theoretical (leading order) distributions are 
shown for all accepted events, while in figures 3b and c these distributions are displayed 
for events with 777,3 + ^4 > 10 GeV and 7773 + 777.4 > 10 GcV respectively. In these three 
cases the difference between experiment and theory can be expressed by the values for 
the chi-squares per degree of freedom, which are 1.8, 0.24, and 4.9 respectively. Visual 
inspection confirms that the difference between experiment and theory is largest for case 
"c" , ma -|- 777,4 > 10 GeV. After adjustment by fit ( the corresponding distributions are 
not shown here ), the chi-squares become 0.9, 0.1 ct 1.3 respectively. Figure 3d is an 
illustration of the theoretical shape differences between qqqq and qqgg events (for all 
accepted events according to ERT), which make our type of analysis possible. 
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4 Discussion 

Figure 2a shows that there are deviations of the experimental distributions from the 
theoretical ones which concentrate in the region of high values of the variable ma + m^. 
Such deviations could also be caused by the absence of higher order effects in the ERT 
simulation. To get an idea of the influence such higher order effects might have, we 
show in figure 2b the results one obtains if one replaces the experimental distributions by 
distributions explicitly containing five-parton events in addition to the 4-parton events 
furnished by ERT. The mixture is adjusted so as to reproduce the observed value of 
1.32 for i?4q. In figure 2c one finds the results for a hybrid calculation, where each of 
the four partons from ERT is followed by a parton shower, before being fragmented into 
particles. These events are also mixed with regular four-parton events so as to reproduce 
the observed value of i?4q. In both cases one expects such calculations to show trends 
associated with higher order effects. In figures 2b and 2c these calculations are represented 
as fictitious data (full circles) to be compared with calculations where such effects are not 
present (histogram). One notices that indeed such higher order effects have a tendency 
to increase with the value of + 7714, but more gradually than the actually observed 
experimental results of figure 2a. Another possibility to explain such results could be the 
appearance of new fermions in the region of 10 to 15 GeV for 777,3 + 7714, i.e. 5 to 7 GeV for 
each of the emitted particles, (see f.ex. refs. [26, 27]). At this stage of the analysis it is 
not possible to discriminate between these possibilities, but it is hoped that this research 
note motivates other groups to verify these results with more precision. 

5 Summary 

The Bengtsson-Zerwas [22] angular correlation, calculated with the second order ERT [3] 
simulation, has been compared with the correlation measured for hadronic four-jet events 
from the disintegration of the Z^. Shape differences between experiment and theory have 
been observed. The differences concentrate in the region of ma + m4 > 10 GeV, where 
ma + m4 is the sum of the intrinsic masses of the two least energetic jets 

References 

[1] H.Jeremie, P.Leblanc and E.Lefebvre, 
Can. Jour. Phys. 84 (2006) 411. 

[2] G.Farrar, Phys.Lett. B 265 (1991) 395. 

[3] R.K.Ellis, D.A. Ross and A.E. Terrano, Nucl. Phys. B 178 (1981) 421. 

[4] T.Sjostrand, Comp. Phys. Comm. 39 (1986) 347; 

T.Sjostrand and M.Bengtsson, Computer Phys. Comm. 43 (1987) 367. 



[5] OPAL collaboration, P.D.Acton et al., Phys. Lett.B 276 (1992) 547. 

[6] J.Andre and T.Sjostrand, Phys. Rev. D 57 (1998) 5767. 

[7] K.Hagiwara and D.Zcppcnfeld, Nucl. Phys. B 313 (1989) 560; 
F.Wackerle, Diplomarbeit Karlsruhe 1994, IEKP-KA/93-19. 

[8] ALEPH collaboration, D.Decamp et al., Phys. Lett. B 284 (1992) 151. 

[9] DELPHI collaboration, P.Abreu et al.,Phys. Lett. B 255 (1991) 466; 
DELPHI coUaboration, P.Abreu et al., Zeit. f. Phys. C 59 (1993) 357. 

[10] DELPHI collaboration, P.Abreu et al., Phys. Lett. B 414 (1997) 401. 

[11] OPAL collaboration, R. Akers et al., Zeit. f. Phys. C 68 (1995) 519. 

[12] OPAL coUaboration, R. Akers et al., Zeit. f. Phys. C 65 (1995) 367. 

[13] ALEPH collaboration, R.Barate et al., Zeit. f. Phys. C 76 (1997) 1. 

[14] OPAL collaboration, G.Abbicndi ct al., Eur. Phys. J. C 20 (2001) 601. 

[15] ALEPH coUaboration, R.Barate et al., Eur. Phys. J. C 27 (2003) 1. 

[16] OPAL coUaboration, K.Ahmet et al., Nucl. Instr. and Meth. A 305 (1991) 275. 

[17] M. Arignon et al., Nucl. Instr. and Meth. A 313 (1992) 103. 

[18] D. Charlton et al., Nucl. Instr. and Meth. A 325 (1993) 129. 

[19] OPAL coUaboration, G.Alexander et al., Zeit. f. Phys. C 52 (1991) 175. 

[20] OPAL coUaboration, K.Ackerstaff et al., Eur. Phys. J. C 2 (1998) 213 

[21] S. Catani, Y. L. Dokshitzer, M.Olsson, G.Turnock and B.R.Webber, Phys. Lett. 
269 (1991) 432. 

[22] M.Bengtsson and P.Zerwas, Phys. Lett. B 208 (1988) 306. 
[23] O.Nachtmann and A.Reiter, Zeit. f. Phys. C 16 (1982) 45. 
[24] H.A. Olson, P. Osland, I. Overbo, Phys. Lett. B 89 (1980) 221. 
[25] JADE coUaboration, S.Bethke et al, Phys. Lett. B 213 (1988) 235. 
[26] E.L.Berger al., Phys. Rev. Lett. 86 (2001) 4231 
[27] E.L.Berger et al, Phys. Rev. D71 (2005) 14007 



